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We investigate the existence of non axisymmetric solutions in the 6-dimensional baby-Skyrme 
brane model. The brane is described by a localized solution to the baby-Skyrme model extending 
in the extra dimensions. Such non symmetric branes have already been constructed in the original 
2 + 1-dimensional baby-Skyrme model in flat space. We generalize this result to the case of grav- 
itating baby-Skyrme and in the context of extradimensions. These non-trivial deformation from 
the axisymmetric shape appear for higher values of the topological charge, so we consider the cases 
of B = 3, 4, where B is the topological charge. We solve the coupled system of the Einstein and 
baby-Skyrme equations by successive over relaxation method. We argue that the result may be a 
possible resolution for the fermion mass hierarchy puzzle. 
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I. INTRODUCTION 



Theories with extradimensions are expected to solve the hierarchy problem and the cosmological constant problem 
£Lj [H • The great challenge with extradimensions is to conciliate the extradimensions with the fact that we observe only 
four dimensions. The typical way is to assume that we live on branes 0. The branes are 4-dimensional subspace of 
the entire extradimcnsional space-time where the matter fields are trapped. One possible extradimensional geometry 
was proposed by Randall- Sundr am (RS), where the extradimensions is warped. The RS brane model is formulated 
in five space-time dimensions [H, Q and shows that the exponential warp factor in the metric can generate a large 
hierarchy of scales. However, such models require a fine tuning between the parameters of the model. Note that in 
the RS model, the brane is infinitely thin. Such a fine tuning is usually not needed if thick branes are considered. 
\ For instance, the brane theories in six dimensions in the abelian strings show a very distinct feature towards the fine- 
tuning and negative tension brane problems Similar compactification was achieved for magnetic monopoles, 
both for positive and negative cosmological constant ■ 

The Skyrme model jTlj is originally a model for the nucleons after a suitable quantization scheme (l2j . There is a 



OO 

1 low-dimensional mimic of the model called baby-Skyrme model. The baby-Skyrme model possesses soliton solutions 
called the baby-skyrmions in two dimensional space [13l4l6i. The warped compactification of the two dimensional 
extra space by such baby-skyrmions was already studied [17( for negative bulk cosmological constant, based on the 
assumption that the cosmological constant inside the three branes is tentatively set to be zero. Addressing the non- 
■ zero cosmological constant inside the branes has been considered for case of the strings [lH and the monopoles 
$H \ Along these directions, we also have studied the baby-skyrmion brane with both positive brane cosmological constant 
and a bulk cosmological constant [2(| . 

The localization of particles such as fermions [13, US HJ and a gauge field [2l| on these branes have also been 
studied. According to the Index theorem a nonzero topological charge implies a zero modes of the Dirac operator 
[22j . The zero crossing modes are found to be the localized fermions on the brane. It follows that the generation of 
the fermions emerges as an effect of the topological charge of the skyrmions with a special quantum number assigned 
by K3. There are different profiles of the zero crossing behavior for different K$, and in our opinion, this might give 
an origin of the finite mass generation fermions. In this paper we begin our analysis for the topological charge B = 3 
because it is the number of particle generations in the fermions. As the straightforward extension, the case of B = 4 
is also examined. 

In our previous studies, we have found an essential requirement for the geometry of the brane solutions: the brane 
should slightly deform from the axisymmetric shape in order to explain the small mass difference between the first 
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two generations of quarks (such as M s — M u ~ 10 2 MeV, which is much smaller than the difference of the first and 
third: M t — M u ~ 10 6 MeV). In [17j , we introduced an effect of the deformation for the matter fields by hand and 
found that about for a slight deformation of the skyrmion the result is in good agreement with the experiment. The 
drawbacks of the analysis are: (i) we need to introduce an arbitrary deformation parameter, (ii) we ignore effect of the 
deformation for the gravity side. Our main concern of this paper is thus to find brane solutions with spontaneously 
broken axisymmetry. The solutions with lower symmetry, such as Za-symmetry without gravity were already found 
using numerical simulation [23| where an energy minimization scheme called the simulated annealing method (SA) (2~H 
was used. For the analysis of the brane solutions, the coupled system of matter field equations and Einstein equations 
has to be treated and the SA scheme is not suitable for the purpose. Therefore, as a first step, we solve baby-Skyrme 
field equations without gravity by employing the standard successive over-relaxation method (SOR) and compare the 
obtained results by the SA as a validation of our code. 

For the brane solutions, we use an ansatz for the line element which is a function of both the two extra coordinates 
(r,9), inspired by the Lewis-Papapetrou ansatz. The resulting coupled partial differential equations are numerically 
solved in terms of the SOR [25j |. 

This paper is organized as follows. In Sec. II, we start with a brief introduction to the baby-Skyrme model and show 
the results using SA. We also present the Eulcr-Lagrangc equations of the non symmetric solutions. The corresponding 
brane model is extensively discussed in Sec. III. In Sec. IV we present the numerical solutions and give a brief summary 
of our results in Sec.V. 



II. THE BABY-SKYRME MODEL AND THE DEFORMED SOLUTIONS 



The lagrangian density of the baby-Skyrme model is given by 

£bab y = y (<9 A/ n) ■ {d M n) - -j(d M n x d N n) 2 - K V(n), (1) 

where n = (774, 712,713) denotes a triplet of real scalar fields with the constraint n • n = 1. The constants ^2,4,0 arc 
parameter with dimension [M], [M _1 ], [M 3 ] respectively, where M is the dimension of mass. The first term in (JTJ is 
nothing but the non-linear a model while the second term is inspired by the quartic term of the 3 + 1 Skyrmc model. 
The potential V is necessary for the stability of the solutions in the sense of Derrick's argument. While the first two 
terms are invariant under a global 0(3) symmetry, n — > Qn,Q 6 0(3), the potential breaks it. Several choices of 
potentials have been investigated. Two 0(2)-symmetric potentials have been intensively studied, namely the old-type 
^(773) = (1 + ra gj fl3ll and the new-type potential ^(773) = (1 — 7i§) [lj|. Potentials with no 0(2) symmetry have 
been studied in 26l. l27j] . In (23|, the authors generalized the old- type potential to ^(773) = (1 + 773 ) s and were able to 
build solutions that deviate from the axial symmetry. The model admits an analytical solution (holomorphic solution) 
for a single soliton for s = 4. The force between two holomorphic solutions is always repulsive so that multi-soliton 
solutions tend to deviate from axial symmetry fH| . 

The topological charge of the baby-skyrmion is defined as 

B = — I ' d 2 x n- (din x d 2 n) (2) 

4w J 

in 2 + 1 flat spacetime. The numerical analysis in [23| has been done by minimizing the following energy functional 
E = J* M jT p dp{^((d p n) 2 + ±(donj 2 )+^±(d p nxd 8 n) 2 + Ko V(n)} (3) 

in the polar coordinate (p, 9) . The method is a kind of the Monte Carlo method in which one generates random 
numbers and properly change the value of the fields n by the numbers so as to fall down the energy. The more 
sophisticated method may be applied to the problem. The simulated annealing method (SA) [24| is the application of 
the Metropolis algorithm which can successfully avoids the unwelcome saddle points. 

The solutions for B ^ 2 exhibit a non-axisymmetric shape for certain parameter ranges of the coupling constant 
«o or of the value of s. We confirm the results of (23| by minimizing the energy per topological charge in terms of 
the SA. In FiglTJ we plot the typical results of the energy densities of the skyrmion with varying kq. As was found 
in [23(1 , at a critical value of Kq, breaking of the rotational symmetry occurs and only the Z2-symmetry survives for 
B = 3. More interestingly, the case of B = 4, the Z2 remains for s = 1, whereas the s = 2 the solution exhibits a 
Z3-symmetry (see FigJ5]). We restrict our numerical analysis for the case of s = 1 except for some results of B = 4. 

In this section, we solve the Euler-Lagrange equations derived from the action ((TJ) and try to reproduce the results 
of the SA. Generally speaking, both results should coincide. However, we need to introduce additional boundary 



FIG. 1: Contour plot of the energy densities of the B = 3 baby-skyrmion for several strength of the parameters \x := kqK2/k\. 
From the left fx = 0.04, 0.06, 0.12, 0.64 with 8 = 1. 



M * '1 



m 

MM 



w <t« 



FIG. 2: Contour plot of the energy densities of the B — 4 baby-skyrmion for the parameter fi — 0.64, s = 1 (left) and 
/i = 0.32, s = 2. 



conditions in order to solve the equations using a SOR algorithm which possibly strongly affect the properties of the 
solutions. 

We use the following ansatz for the scalar triplet n: 



n = (sin F(p, 9) cos Q(p, 9), sin F(p, 9) sin Q(p, 6), cos F(p, 9)), 

where Q(p, 9) := Oo(p, 0) + ra0 where Qo(p, 9): S 1 h-> 5 1 is homotopic to the constant map. 
The resulting Euler-Lagrange equations are of the form 



(4) 



r 2 d 2 F + d 2 F + rd r F - - sin2F(r 2 (d r Q) 2 + (d e O) 2 ) 



sin 2 (d r Fd e Q - d e Fd r Qjd e e + d 2 F (d e Q) 2 + d 2 g F (d r Q) 2 + (d r Fd e B + d 9 Fd r Gjd r d e e 

■ + ^sm2Ffd r Fd 9 e-deFdrO) - ^ = , (5) 



r 2 5 2 9 + d 2 e Q + rd r Q + cot F(r 2 d r Fd r Q + d e Fd e &) 

(drFdeQ - d e Fd r QjdeF - d 2 .Q {d F) 2 - d 2 e Q (d r F) 2 - 



d r Fd e e + d e Fd r Q )d r d e F 



d 2 Fd e Fd e Q + d 2 Fd r Fd r Q + 2d r d e &d e Fd r F 



(6) 



where we defined dimensionless quantities according to r := y ^2/ K4P, M : = KQK4/K 2 ,. 




FIG. 3: Cartoon figure of the arrow plot of (ni, n-i) in the case of three centered solution of the B — 3. 
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FIG. 4: The energy per topological charge of the B — 3 as a function of the "delocalization " ro, in the case of fj. = 0.08. From 
the result, we extract ro = 2.6 as the true distance of the centers. 



The static energy ^ reduces to 

^static = E 2 + E4 + Eq, 



E 2 := j d9 J rdrj U(d r F) 2 + sin 2 F(d r e) 2 ^j + ^ ([d e F) 2 + sin 2 F(dgO)'- 
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and the topological charge JI) reduces to 



d6 j rdr— sin 2 F ( d r Fd e O - d e Fd r Q 



d6 I rdr b(r,0); b(r,9) := -- — [d r (cos Fd g Q) - d e (cos Fd r Q)]. 
We introduce a deformation factor A in order to quantify the deformation of the solutions, defined as 

n 2\ 



a 2 27r 

A = — 
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rdr 6(r, 6) 
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FIG. 5: The profiles F(r, 9), 0(r, 8) of the B = 3 as solutions of the Euler-Lagrange equations with /i = OJ 
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FIG. 6: The energy density of the B = 3 solution. 
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The B = 3 solution 



The Z2-symmetry imposes that the equations are invariant under the transformations 

6»^6» + 7r: F(r, 61) -> F(r, 9 + n), Q (r, 9) -> 9 (r, + tt) 

and at the axis of symmetry = 0, 

tf^-tf: F(r,fl)->F(r,-ff), e o (r,0) -> -0 o (r,-0). 

Therefore we restrict our numerical calculations to the half-plane defined by ^ r ^ oo and — £ ^ ^ 5- . 
We tried using the following boundary conditions, 

F(O,0)=O, 5 r e(r,6»)| r=0 = 0, F(oo,9)=ir, d r Q(r, 9)\ r=OQ = 0, for 

and 

¥M)|«= ±f =0, e(r,±|) = ±|7r > for O^r^oo. 

Although they are sufficient to get the solutions with axisymmetry (i.e., the left panel of FigQJ, they are not suited 
for the symmetry breaking solutions and additional conditions seem to be required to reproduce the results. The 
reason is that for a critical value of /x, peaks appear and grow at a finite distance from the origin (see Fig[T]). These 
points are singular points where the direction of the component arrow of n is frustrated. The geometric structure is 
schematically drawn in Figj3] one easily see that for r < ro (here r is the distance from the origin), the geometry looks 
like the one center hedgehog, while for r > r the arrow winds three times around the origin. We take this property 
of the solution geometry into account by using the following boundary conditions (|12[) . 



(10) 



(11) 



(12) 



(13) 



d e F(r,0)\9=±s 
f(ro,±f) = 



0, for ^ r < r , r < r ^ oo, 



(14) 
(15) 
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FIG. 7: The figures of the arrow plot of n of the B = 3 solutions in the case of the one centered, ro = solution (left) and 
the three centered, ro ~ 3.0 solution (right) for /i = 0.08. The length of the arrow is relative to the norm of the vector, i.e. 



and 

f ±| ^ r < r 

e ( , ± f)= do, 

^ ig 71 " r ° — r — 00 

where ro is a constant describing the center distance of the solution constituent solitons. Note that when ro — > 
the boundary conditions (fl4)) - (fT6|) coincide with the case of the axisymmetry (fT3|) . One can easily check that these 
modified boundary conditions correctly reproduce B = 3 (see appendix). 

The equations with the appropriate boundary conditions can be solved numerically by the SOR method. In order 
to get the non-axisymmetric solution, we employ the following procedures. We begin the SOR for a fixed ro with 
following initial conditions which satisfies the boundary behavior of the axisymmctric solutions 

F(r,fl)|i niU Bi = 7r(l-e- ArS ) ) 0(r, 0)| initial = 30. (17) 

For changing ro, we repeat the computation and finally reach the true value of ro by minimizing the energy par 
topological charge E/B, which realizes the true boundary condition for the system. In terms of the procedure, we 
successfully find the non-axisymmetric shape of the solutions. 

In FigU]we plot a typical E/B as a function of ro. In this case, we could choose ro = 2.6 as a minimized energy 
per unit charge. Fig f5] shows a numerical result of the profile functions F(p,0),Q(p,9). In FigJHl we present the 
corresponding topological charge density for several potential parameter \i. One easily sees that for increasing values of 
fi, the peaks grow and a new peak appears at the origin. Note that the F, O profiles contain singularities/discontinuities 
corresponding to the centers of the constituents. In Fig|7]we plot the configurations of n, i.e., the vectors (ni, riz) on 
the plane for the case of ro = 0.0 and ro ~ 3.0 which corresponds to minima of the case of n = 0.08. One easily see 
that these are perfectly regular. 

In Fig|5]we plot the behavior of the deformation parameter A for changing the strength of the potential. For larger 
/Lt, the deformation monotonically grows. Note that there is a critical point of the deformation of the solution: for 
fi ~ 0.04, deformed solutions suddenly appear. The error-bar is due to the fact that the results involve a systematic 
uncertainly for the determination of ro . 

The solution of the Euler-Lagrange equation should coincide with result of the energy-minimization scheme. In 
FiglHl we compare the energy E and also the energy per topological charge E/B in terms of the SOR and of the SA 
for the potential parameters /i. Both methods contain uncertainties in several reasons, in particular the SA always 
tends to underestimate the topological charge as well as the energy. As a result, the energy of the SA is lower value 
than that of the SOR while the energy per topological charge is slightly higher, but the results are in good agreement 
within ~ 5%. 

For higher charge, solutions breaking the axisymmetry appear even for very small fi. 
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FIG. 8: The deformation parameter A of the B = 3 solution as a function of ji. 




FIG. 9: The comparison of the static energies (left) and the energies per the topological charge (right) of the B = 3 solutions 
for the successive over-relaxation (straight line) and the simulated annealing (dotted line). 



B. The B = 4 solution 



The cases with a Z2-symmetry are treated essentially in the same way as in the B = 3 case. Therefore, here we 
concentrate on the case of Z3-symmetry. A typical result of the energy density is shown in Fig. [2] 
For the Z3-symmetric case, the equations must be invariant under the transformations 

2n 2-7T 2n 

9^6+ — : F(r,6)^F(r,6+—), Q (r, 9) -> o (r, 9 + — ) (18) 

and if we assume that one of the peak is located on the x-axis (9 = n), the equations are invariant under the 
transformations 

8^-9: F(r,6)^F(r,-9), 9 (r, 9) -> -6 (r, -9). (19) 

Therefore we perform our numerical calculations on the one third- plane defined as ^ r ^ oo and — ^ ^ 6 ^ -| . The 
functions on the whole plane can be estimated in terms of the above symmetry. 



FIG. 10: The profiles F(r,8), Q(r,6) of the B = 4 as solutions of the Euler-Lagrange equations with s = 2, fj, = 0.01. 



We use the following boundary conditions: 



and 



F(0,6) = 0, d r G(r,6)\ r = = 0, F(oo,0)=tt, d r G(r,6)\ r=00 = 0, for (20) 

%F(r,0)| e=±¥ =0, e(r,±£)=±y, for O^r^oo (21) 

for the axial symmetric solutions and when the constituent soliton centers are located on the line of 9 = 7r , we 
take the singularities of the centers into account by imposing 

F(O,0)=O, a r e(r^)| r=o = 0, F(oo,0)=7r, 3 r e(r, 0)| r=oo = 0, for (22) 

a e F(r, 0)| fl=±i =0, for ^ r < r , r < r ^ oo, (23) 

^(r o ,±J) = (24) 



and 



±| ^ r < r Q 

>±l)={ A ( 25 ) 

±— roSrS oo. 
o 

Again for ro — > the boundary conditions reduces to the axisymmetric case. 
Fig [TU1 shows a numerical result of the profile functions F(p, 9), 0(p, 9). 

III. THE GRAVITATING BABY-SKYRME MODEL IN SIX DIMENSIONS 

We are now are ready to include gravity to the baby-Skyrme brane model. In this paper, we will focus on the B = 3 
but the method is straightforwardly extendable to the cases of different B. We also show the result for B = 4 as an 
example. 

The total action for the gravitating baby-Skyrme model in six dimensions is of the form S = S gra:v + Sbaby with 

Sgrav = j <?xV=-g {^-R - A(6)J , (26) 

Sbaby = / (FxyZ-gCbaby, (27) 
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where A( 6 j is the bulk cosmological constant, and X(6) — 87rG( 6 ) = 87r/M^. 

The coefficients ^2,4.0 hi Eq. (|27[) are the dimensionful coupling constants in the gravitating baby-Skyrme model 
with dimensions [A (6) ] = M 6 , [x (6) ] = M~ 4 , [« 2 ] = M 4 , [k 4 ] = M 2 [«o] = M 6 . 

A. The Ansatz 

We consider the case of axial symmetry breaking in the extra dimensions, and we use the following metric ansatz 
inspired by the Lewis-Papapetrou ansatz in 4-dimensions 



ds 2 = M\p, 6)g$dx»dx» - ^Jj^ + P 2 dB 2 ) 



(28) 



where p £ [0, 00) and 9 G [— 7r, k] are the coordinates associated with the extra dimensions. We further model a 
cosmological constant on the brane by considering the following form of the four dimensional subspace (described by 

J 4 1 in Pn " 

g^v m &q. 



ds 2 {i) = gftdx^dx" = di 2 - <%e 2ff(t WcfaJ'. 
H(t) is a function of the time coordinate t which describes the inflation in the four dimensions. 



(29) 



B. Field equations of the model 

In order to rewrite the system in terms of dimensionless quantities, wc dehne the dimcnsionless variable 



V «4 V K 2 



(30) 



and dimcnsionless parameters 



X(6)«2, P:=A (6) — , p:=n — . 

K 2 K 2 



(31) 



After straightforward computations, the Einstein equations reduce to 
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M 2 



a(re - /?), 



(32) 



(33) 



(34) 



(35) 



where we used the notation T~L,t '■= dtH{t) and H.t,t '■= d 2 l~L(t). Let us note that (|32l) and (|33|) are the four-dimensional 
components of the Einstein equation, while (|34|) and (|35|) are the extra-dimensional components. The components of 
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■((d r F) 2 + sin 2 F{d r Q) 2 ) 



M2 ((d e F) 2 + sin 2 F(dgQ) 2 ) - ^ 



the dimcnsionless energy- momentum (EM) tensor in Eqs. (|32l) - (|35p are given by 
M 2 

— ((f)„ F\ 2 -i- sin 2 F<7Lft W — ((r)n F') 2 4- sin 2 F(f)nP\\ 2 \ — sir 2 

2L 2 

M 2 

Tr ~ 2L 2 ^ ' * K " r "' ' 2r 2 L 2 ^^ ' ' ~~ ~ y "°~' ' ' 2r 2 L 4 

M 2 M 2 M 4 

ro = ~^((d r F) 2 + sm 2 F(d r e) 2 ) + ^ 7I ((a e F) 2 + sin 2 F(a e) 2 ) + ^ 7I sin 2 F(a [rJ Fa e] e) 2 - M y(n), 



2r 2 L 2 
M 2 



((d e F) 2 +sm 2 F(d e Qf 



2r 2 L 4 
M 4 



sin 2 F(d [r Fd e] ey - nV(n), 
sin 2 F(d [r Fd e] Qf - //y(n), (36) 



2r 2 L 2 



2r 2 L 4 



where d^Fdg^Q. := <9 r .F<9 e e - d e Fd r O. 

From (|32|) . (|33l) one directly sees that the function "H(i) must be linear in time: 

«(f) = Hot (37) 

where Ho is a constant called Hubble parameter. For later convenience, we introduce a dimcnsionless parameter 

(38) 



K2 



We finally get the following equations by using a suitable linear combination of (|32p - ([55l) 
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For the matter fields, the equations are given by 
M 2 r „, „ / . d r M 1 
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where d{ r FdgyQ := d r Fdg<d + dgFd r Q and d r g := d r dg. 

A useful quantity for the gravity-side is the rescaled Ricci scalar which is given by 



R(r,6) 



M 2 r d 2 M 

G^rr- + 14 
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(39) 



(40) 



(41) 



(42) 



(43) 



C. Boundary conditions 



For the Z2-symmetry, the Einstein equations ([32j) -(|36 |h and the matter field equations (|4T j) .(|42 |) must be invariant 
under the transformations 

e^e + ix: F(r,6) -+ F(r,6 + 7r), 9 (r, 6) 9 (r, 6 + tt), M(r, 0) -H- M(r, + vr), L(r, 6) -H- L(r, + vr) (44) 
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and when we set 9 = as the axis of symmetry, the equations mast be invariant under the transformations 

9^-9: F(r,6) -> F(r,-6), Q (r,6) -Q (r,-6), M(r,9) -H- M(r,-9), L(r,9) -H- L(r,-6). (45) 

Thus as in the case of non gravitation, it is sufficient to restrict the integration domain to the half-plane defined by 
^ r ^ oo and — ^ ^ 9 ^ ^. For the matter fields, we employ the boundary conditions (fT12j) and (p^|) - (PTB]) . 
For the metric fields, regularity at the core of the brane leads to the following conditions 



Af(O,0) = l d r M(r,9)\ r=0 = 0- L(0,9) = l, d r L(r 7 0)\ r=o = 0, for --^9^-. 



(46) 
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Thus it is likely that the solutions satisfy the following additional boundary conditions 

d e M(r,6)\ e=±f = 0, d 9 L(r,9)\ e=±f = 0, for ^ r <: oo. (47) 

Note that we do not need to impose any conditions for the metrics about ro because the local spatial structure of the 
brane is described by only the matter fields. 

D. The potential 

As we already mentioned, in (23j the authors have introduced a two parameters class of potential 

V = /i(l + n 3 ) s < s < 4. (48) 

In the original approach, the authors used dimensionful parameters (which is essentially our K4/K2) denoted by n 2 
and minimized the energy for changing values of s, k 2 with fixed \i. The main difference in our analysis is that we 
vary the parameters of the potential, i.e. /1 and successively solve the equations instead of minimizing the energy. 

IV. NUMERICAL RESULTS 

In [20j ] . we have done the analysis for the same model but have assumed axisymmetry from the start. In this case 
the problem reduces to a system of coupled ordinary differential equations. In the present case we have to treat the 
four coupled partial differential equations (j3"2")) -([4*2" j) with the boundary conditions p^)) . p^| - p3|) by using the SOR 
method. 

We use the following form 

F(r,0)| initia] := tt(1 - e^ 3 ), 0(r, 0)| initial := W 

M(r, ^1™^ = L(r,0)| initia i = 1 (49) 

as the initial profile of the relaxation scheme. 

Note that the conditions we use, together with the initial profile (|49|) leads to axisymmetric solution. However, 
when the deformation is relatively small, they are approximativcly consistent with the boundary conditions of the 
non-axisymmetric case. Thus, we begin with ro = and find the solution. By using this result as a initial profile, 
the computation for the small finite ro is performed. We repeat the procedure until we get the solutions for sufficient 
large value of ro. For the continuous change of ro, we can find minimum of the energy per topological charge. Since 
the true solution should exhibits the minimum energy per charge, we employ the ro corresponding to the minimum 
energy per charge for the boundary conditions. 

Since the boundary conditions for the metrics are imposed only at the origin, the Einstein equations essentially 
are the initial value problem. Usually it is not straightforward to adapt the SOR to the problem. Thus we employ 
the following procedure. We put an initial guess of the form of L, M at the infinity and solve the equations for that 
condition. If the condition is inappropriate, the solution surely exhibits the non-trivial singularity at some location, 
then we explore the proper boundary condition until such singular behavior disappears. 

In FigfTTJwc show a solution for F, 9. The metric functions L, M are shown in Fig[T2]for both sign of the cosmological 
constant /3. We plot the scalar curvature corresponding to the solutions in FigfT51 In both cases, there are two distinct 
peaks together with the remaining slight ridge. The result for the fi = 0.08; if we choose larger value of /i, the peaks 
grow and finally end up in three independent peaks corresponding to the energy density of the baby-skyrmion. 

The Ricci scalar for a typical solution with the Z3 symmetry is presented in Fig|TJ] We clearly see three distinct 
peaks and for larger value of //, the fourth peak grows at the origin until we finally get a four-centered solution. 

V. SUMMARY 

In this paper, we have obtained the non-axisymmetric brane solutions in the B = 3,4 baby-Skyrme model with 
the old-baby potential. We treated the problem in two steps: we first computed the solution using the SA in order 
to gain intuition on the boundary conditions. Then we solved the four coupled partial differential equations by using 
the standard SOR by using suitable boundary conditions. The axial symmetry breaks to Z2 symmetry and the 
deformation is a one-parameter family of the strength of the potential fi. 
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= 0.05 P = -0.05 

FIG. 13: The Ricci scalar curvature of the B = 3 solutions for a = 0.005, 7 = 0.0 . 




FIG. 14: The Ricci scalar curvature of the B = A for a = 0.005, P = 0.0, fi = 0.01, s = 2. 

Recently the authors have found at a related model with the present study that there are holomorphic solutions for 
special choice of the potentials [28|, H^] . The potentials are a class of old-baby potential for B = 1 while for all other 
B they become a class of new-baby. Furthermore, the axisymmetry is a good symmetry for the new-baby potentials 
while for the old-baby the symmetry is inevitably broken. The most crucial thing for the deformation is thus whether 
the value of the potential at the origin is finite or not. 

Study of the linear stability of the present model has been done in [20I ] within the assumption of the axial symmetry. 
Apart from the axial symmetry, we also need to take into account the fluctuation of 0. It seems straightforward but 
surely cumbersome. For the level of the classical solution, our solution is stable for the variation of the size moduli 
(see Fig|4]), so we expect that our solution will be stable against any general perturbation as well. 

Of course there are several possibilities of the potential for obtaining the anisotropic solutions. For example, in 
[27l ] the authors introduced a more sophisticated potential V = |1 — (ni + iri2) N \ 2 (l — n 3 ) (with the integer N ^ 2). 
However the potential is essentially of old-type and the solution tends to split the number of B constituents. An 
interesting feature of such models is the fine structure: the zero of the potential are at the infinity (F = 0) as well as 
the points F = £ and = ^ (n:integer). As a result, the solution split B component where each has NB peaks. 

An intcrsting question raised by this work is to see the effects of the the spatial structure of the branes on the 
mass of the brane localized fermions (quarks, leptons). The fermions coupled with the baby-skyrmions with Yukawa- 
coupling exhibit non-trivial doubly degenerated plus isolated states which could be identified as three generations of 
the quarks/leptons. The model has weak isospin symmetry and the isodoublet of the quarks are degenerate. We 
expect that the brane slightly deforms from the axisymmetry and reproduce the slight splitting the mass of the first 
and second levels [13] • Now we got the intuition that the level splitting of the generations is owed by the coupling 
constant of the old-baby potential \i. More thorough analysis of the property of the localized fermions in these funny 
structures of the brane will be certainly very interesting. 
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Appendix A: Topological charge 

The topological charge defined by ([2]) is easily computed for the boundary conditions (fl~2|) . (fT4^ - (fl~6j) . From 

B = -— / d9 dr [d r (cosFd g Q)-dg{cosFd r Q)] 
2tt J-iL Jo 

= — — — / dB \cosFd e e}°° + — / dr [cosF^el 2 „ 
1 rf 

— / d6 (cos F\ r=oo d0<d\r=oo - cosF| r=0 9e9| r= o) 



2- 
1 

'2n Jo 



dr (cosFU = a.9 r 6|e = 5 - cosF| e= _-9 r e| e= _: 



= -^((-!) • 37T - 1 • TT) + ^(1 • 7T - 1 • (-7T)) 

= 3 (Al) 
where we used the relation ±-^d r Q\g = ±z. = S(r — r ). 
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